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We compare nuclear and neutron matter predictions based on two different ab initio approaches to 
nuclear forces and the nuclear many-body problem. The first consists of a realistic meson-theoretic 
nucleon-nucleon potential together with the relativistic counterpart of the Brueckner-Hartree-Fock 
theory of nuclear matter. The second is based on chiral effective field theory, with density-dependent 
interactions derived from leading order chiral three-nucleon forces. We find the results to be very 
close and conclude that both approaches contain important features governing the physics of nuclear 
and neutron matter. 



I. INTRODUCTION 



Nuclear matter is a convenient laboratory to test nuclear forces and many body theories. In particular, the 
equation of state (EoS) of extremely neutron-rich matter has attracted considerable attention lately because of its 
broad applications, ranging from the structure of rare isotopes to the properties of neutron stars. 

Constraints on the properties of nuclear matter, symmetric or isospin-asymmetric, can be obtained from a va- 
riety of experiments, such as measurements of nuclear binding energies (including isobaric analog state energies), 
parity-violating electron scattering, neutron skin thickness measurements, nucleus-nucleus collisions, and astrophysi- 
cal observations. For a recent review of available constraints, particularly on the nuclear symmetry energy, the reader 
is referred to Ref. [T]. 

Diverse theoretical frameworks have been employed to describe the properties of nuclear and neutron mat- 
ter. They include: phenomenological approaches, both relativistic and non-relativistic; the Brueckner-Hartree-Fock 
(BHF) method, typically implemented with three-nucleon forces (3NF); variational approaches; the relativistic Dirac- 
Brueckner-Hartree-Fock (DBHF) approach; and chiral effective field theories. Predictions are model-dependent, par- 
ticularly at the higher densities, where constraints are scarse and less stringent. 

In our previous work with both symmetric nuclear matter (SNM) and isospin-asymmetric nuclear matter (IANM), 
we have relied on the DBHF approach, together with a relativistic meson-theoretic nucleon-nucleon (NN) potential 
which uses the pseudovector coupling for the pion [2] . The latter choice was motivated by the considerations we outline 
next. Already when QCD (and its approximate symmetries) were unknown, it was observed that the contribution 
from the nucleon-antinucleon pair diagram, Fig. [I] becomes unreasonably large if the pseudoscalar (ps) coupling is 
used for the pion, leading to very large pion-nucleon scattering lengths [H 0] . We recall that the Lagrangian density 
for pseudoscalar coupling of the nucleon field (tjj) with a pseudoscalar meson field (cf>) is 

C ps = -ig ps tpj 5 ip(p. (1) 

On the other hand, the same contribution shown in Fig.JT] is heavily suppressed when the pseudovector (pv) coupling 
is used instead (a mechanism which became known as "pair suppression" ) . The reason for the suppression is that the 
covariant derivative in the pseudovector Lagrangian, 

£ P v = -^^757 M W , (2) 

nips 

generates a vertex that is proportional to momentum (leading to a weak coupling for low momenta) and, thus, explains 
the small value of the pion-nucleon scattering length at threshold. Non-linear realizations of chiral symmetry [3] can 
further motivate a preference for the pseudovector coupling. 

Back to the many-body problem, the main strength of the DBHF approach is its inherent ability to account for 
important three-body forces through its density dependence. In Fig. [2] we show a three-body force originating from 
virtual excitation of a nucleon-antinucleon pair, known as "Z-graph" . Notice that the observations from the previous 
paragraph ensure that this diagram, as well as the corresponding diagram at the two-body level, Fig. [T] is moderate 
in size when the pv coupling, Eq. ([2]), is used. (Hence, the importance of using pv potentials in relativistic nuclear 
structure calculations.) The main feature of the DBHF method turns out to be closely related to the 3NF depicted 
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FIG. 1: Contribution to the NN interaction from virtual pair excitation. Upward- and downward-pointing arrows represent 
nucleons and antinucleons, respectively. Dashed lines denote mesons. 



FIG. 2: Three-body force due to virtual pair excitation. Notation as in the previous figure. 



in Fig. [2] as we will argue next. In the DBHF approach, one describes the positive energy solutions of the Dirac 
equation in the medium as 



u*(p,A) = 



El 



2m* 
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1 

(T-f) 
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XX, 



(3) 



where the nucleon effective mass, m* , is defined as m* = m + Us, with Us an attractive scalar potential. (This will 
be derived belowj It can be shown that both the description of a single-nucleon via Eq. ^ and the evaluation of the 
Z-diagram, Fig. ^1 generate a repulsive effect on the energy/particle in symmetric nuclear matter which depends on 
the density approximately as 



AE 



8/3 



(4) 



and provides the saturating mechanism missing from conventional Brucckncr calculations. 

The approximate equivalence of the effective-mass description of Dirac states and the contribution from the Z- 
diagram has a simple intuitive explanation in the observation that Eq. |3]), like any other solution of the Dirac 
equation, can be written as a superposition of positive and negative energy solutions using free nucleon masses. On 
the other hand, the "nucleon" in the middle of the Z-diagram, Fig. [2] is precisely a superposition of positive and 
negative energy states. In summary, the DBHF method effectively takes into account a particular class of 3NF, which 
are crucial for nuclear matter saturation. 

The conventional BHF theory, together with meson-exchange 3NF, can also describe the saturation of nuclear 
matter in a satisfactory manner, although consistency between the parameters of the two- and the three-body forces 
can be problematic [5]- 

Ideally, one wishes to derive nuclear forces from the fundamental theory of strong interactions, QCD. Such task 
is however not feasible, due to the non-perturbative nature of the theory in the low-energy regime typical of nuclear 
physics. An alternative way is to respect the symmetries of the QCD Lagrangian while retaining the traditional 
degrees of freedom relevant to nuclear physics, namely nucleons and pions [5J. This is the philosophy of chiral effective 
field theorie (EFT), which has become popular in recent years. 

Effective held theories allow for a systematic expansion in powers of the momentum (or the pion mass) known 
as chiral perurbation theory (ChPT) such that, at any order, the irreducible two- and many-body diagrams to be 
included are precisely defined. An extensive review on the subject, including a comprehensive list or references, can 
be found in Ref. [7] . 
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Chiral EFT has validity only up to the chiral symmetry breaking scale of A x ss 1 GeV. Thus, the low-momentum 
expansion has a limited range of applicability. This is another reason why relativistic meson theory has been considered 
more suitable for applications to dense systems, where high momenta are involved due to the high Fermi momenta. 
To remedy this problem, ways to extend EFT predictions to higher densities are employed, such as parametrizing the 
available predictions with accurate fitting functions which are then used to predict the properties of dense astrophysical 
systems. 

It is the purpose of this paper to conduct a comparison between nuclear and neutron matter (NM) predictions 
obtained with two ab-initio approaches: 

1. The DBHF method as briefly outlined above, using the Bonn B meson-exchange potential [8 . 

2. A quantitative chiral NN potential, for which we choose the Idaho N 3 LO potential of Ref. [3] together with 
chiral three-nucleon forces at N 2 LO. Details are given in the next section. 

We would like to explore how the DBHF phenomenology compares with the predictions of chiral EFT at N 3 LO 
plus leading 3NF, in both symmetric and neutron matter. 

We will start with a comparison at the two-body force level. The meson-exchange side of this part of the comparison 
will be represented by the conventional BHF model. In the meson-exchange sector, realistic saturation properties will 
be generated through the DBHF model, whereas, on the chiral side, the contribution from leading 3NF will be 
included. We wish to compare the size and density dependence of the saturating effects in each scheme. We will 
compare predictions of the equation of state and the symmetry energy in both cases. Neutron star masses and radii 
will also be addressed. 

This paper is organized as follows: In Section II, we review the main points of the DBHF calculation leading to 
the energy per particle in nuclear matter with arbitrary degree of isospin asymmetry. A brief outline of chiral two- 
and three-nucleon forces is contained in Section III. Our findings and conclusions are presented in Sections IV and V, 
respectively. 



II. RELATIVISTIC MESON-EXCHANGE POTENTIAL AND DBHF APPROACH 



As stated in the Introduction, the starting point of our many-body calculation is a realistic NN interaction which 
is then applied in the nuclear medium without any additional free parameters. 

Relativistic meson theory is an appropriate framework to deal with the high momenta encountered in dense matter. 
In particular, the one-boson-exchange (OBE) model has proven very successful in describing NN data in free space up 
to high energy and has a good theoretical foundation. The OBE potential is defined as a sum of one-particle-exchange 
amplitudes of certain bosons with given mass and coupling. In general, six non-strange bosons with masses below 1 
GeV/c 2 are used. Thus, 



E 



v, 

CZ=7T ,T},p,UJ,5,<7 



(5) 



with 7r and 77 pseudoscalar, a and 8 scalar, and p and uj vector particles. For more details, see Ref. [8]. 

Among the many available OBE potentials, some being part of the "high-precision generation" [TDJ [TT], we seek a 
momentum-space potential developed within a relativistic scattering equation, such as the one obtained through the 
Thompson [T2] three-dimensional reduction of the Bethe-Salpeter equation [T3] . 

Having summarized in the Introduction the main DBHF philosophy, for completeness we now proceed to review 
the main aspects of our approach and the various approximations we perform through the application of the DBHF 
procedure. The equations we present are those suitable for isospin-asymmetric nuclear matter (IANM), since they 
naturally contain both the cases of SNM and NM. 

We start from the Thompson [12] relativistic three-dimensional reduction of the Bethe-Salpeter equation [13] . The 
Thompson equation is applied to nuclear matter in strict analogy to free-space scattering and reads, in the nuclear 
matter rest frame, 

gij(cf, q, p, (4')o) = <j(<f , q) 

d3R ... r* &\ m i m i Qio(K,P) 



*J WF*' (e?J„ - i{ j (P,K) g ' 1 ' iK '*' P ' ( '* J> ° ) ' <« 

where gij is the in-medium reaction matrix (ij=nn, pp, or np), and the asterix signifies that medium effects are 
applied to those quantities. Thus the NN potential, v*^ is constructed in terms of effective Dirac states (in- medium 
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spinors) as explained above. In Eq. q, <f , and K are the initial, final, and intermediate relative momenta, and 
E* = J (m*) 2 + K 2 . The momenta of the two interacting particles in the nuclear matter rest frame have been 
expressed in terms of their relative momentum and the center-of-mass momentum, P, through 



and 



The energy of the two-particle system is 



P = h+k 2 (7) 



K=^. (8) 



i (P,K) = e*(P,K) + e*JP,K) (9) 



and (e^o is the starting energy. The single-particle energy e* includes kinetic energy and potential energy contri- 
butions (see Eq. (23) below). The Pauli operator, Qy, prevents scattering to occupied nn, pp, or np states. To 



eliminate the angular dependence from the kernel of Eq. ^ , it is customary to replace the exact Pauli operator with 
its angle-average. Detailed expressions for the Pauli operator and the average center-of-mass momentum in the case 
of two different Fermi seas can be found in Ref. [T4"] . 



With the definitions 



G ^7^^lS^ (10) 



and 



V ■ = — v ■ — (11) 



one can rewrite Eq. (|6| as 

G l3 {^q,P, (e* j ) ) = Vf j (q',q) 

^^■ g} ^-S,g) G ''^ P ^ M ' (12) 

which is formally identical to its non-relativistic counterpart. 

The goal is to determine self-consistently the nuclear matter single-particle potential which, in IANM, will be 
different for neutrons and protons. To facilitate the description of the procedure, we will use a schematic notation for 
the neutron/proton potential. We write, for neutrons, 

U n = U n p -\- U nn , (1^) 

and for protons 

Up = U pn + Upp , (14) 



where each of the four pieces on the right- hand-side of Eqs. ( l~3p4 ) signifies an integral of the appropriate G-matrix 



elements (nn, pp, or np) obtained from Eq. (12). Clearly, the two equations above are coupled through the np 
component and so they must be solved simultaneously. Furthermore, the G-matrix equation and Eqs. ( l~3|[l4 | are 



coupled through the single-particle energy (which includes the single-particle potential, itself defined in terms of the 
G-matrix). So we have a coupled system to be solved self-consistently. 

Before proceeding with the self-consistency, one needs an ansatz for the single-particle potential. The latter is 
suggested by the most general structure of the nucleon self-energy operator consistent with all symmetry requirements. 
That is: 

UM = UsM + loK tl ( P ) - 7 • pUvAp) . (15) 

where U$\i and Uy,i are an attractive scalar field and a repulsive vector field, respectively, with Uy i the timelike 
component of the vector field. These fields are in general density and momentum dependent. We take 

Ui(p)~U s ,i(p) + ~foU v Ap), (16) 
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which amounts to assuming that the spacelike component of the vector field is much smaller than both Us,i and U v 



v,i- 



Furthermore, neglecting the momentum dependence of the scalar and vector fields and inserting Eq. ( 16 ) in the Dirac 
equation for neutrons/protons propagating in nuclear matter, 



{ipP^ - mi- Ui(p))ui(p, A) = , 
naturally leads to rewriting the Dirac equation in the form 

(7/i (pT - m*)ui(p, A) = , 
with positive energy solutions as in Eq. ([3]), m* = m + Usj, and 

(p°r=p°-u v j P ). 



(17) 
(18) 
(19) 



The subscript signifies that these parameters are different for protons and neutrons. 

As in the symmetric matter case |15j . evaluating the ex pectatio n value of Eq. (16) leads to a parametrization of 
the single particle potential for protons and neutrons (Eqs.( l3|l4 |) in terms of the constants Us.i and U Vi which is 
given by 



Also, 



u i (p) = r ^<p\u i (p)W>= r ^u s , 



u i(p) = E E G iM$) 



V.i ■ 



(20) 



(21) 



The kinetic contribution to the single-particle energy is 



which, along with Eq. (|20|), allows the self-consistent determination of the single-particle potential as explained below. 

(22) 



T AP) = < PI7 ' P + m \P > 



E* 



- P 2 



and the single-particle energy is 



The constants m* and 



e*(p)=T i (p) + U i (p)=E* 



U ,i = U s , t + Uy 



V.i 



(23) 
(24) 



are convenient to work with as they facilitate the connection with the usual non-relativistic framework [16j . 

Starting from some initial values of m* and Ua,i, the G-matrix equation is solved and a first approximation for 
Ui(p) is obtained by integrating the G-matrix over the appropriate Fermi sea, see Eq. ( |21[ ). This solution is again 
parametrized in terms of a new set of constants, determined by fitting the parametrized t/j, Eq. (20), to its values 
calculated at two momenta, a procedure known as the "reference spectrum approximation" . The iterative procedure 
is repeated until satisfactory convergence is reached. 

Finally, the energy per neutron or proton in nuclear matter is calculated from the average values of the kinetic and 
potential energies as 



d = — < % > <Ui> -m. 
The EoS, or energy per nucleon as a function of density, is then written as 

_/ \ Pn&n T" PpCp 
e{Pn,Pp) = , 



or 



e(fcp, a) 



(1 + a)e n + (1 - a)e p 



(25) 



(26) 



(27) 



Clearly, symmetric nuclear matter is obtained as a by-product of the calculation described above by setting a—0, 
whereas a=l corresponds to pure neutron matter. 
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III. CHIRAL TWO- AND THREE-NUCLEON FORCES 

Ideally, one wishes to base a derivation of the nuclear force on QCD. However, the well-known problem with QCD 
is that it is non-perturbative in the low-energy regime characteristic for nuclear physics. For many years this fact was 
perceived as a great obstacle to a derivation of nuclear forces from QCD-impossible to overcome except with lattice 
QCD. The effective field theory concept has shown the way out of this dilemma. One has to realize that the scenario 
of low-energy QCD is characterized by pions and nucleons interacting via a force governed by spontaneously broken 
approximate chiral symmetry. This chiral EFT allows for a systematic low-momentum expansion known as chiral 
perturbation theory (ChPT). Contributions are analyzed in terms of powers of small external momenta over the large 
scale, (Q/A x ) u , where Q is generic for an external momentum (nucleon three-momentum or pion four-momentum) or 
pion mass, and A x s» 1 GeV is the chiral symmetry breaking scale ('hard scale'). (See Ref. j7] and references therein.) 

The past fifteen years have seen great progress in applying ChPT to nuclear forces. As a result, NN potentials 
of high precision have been constructed, which are based on ChPT carried to next-to-next-to-next-to-leading order 
(N 3 LO). We will apply here the chiral NN potential of Ref. [5] which uses a cutoff, A, equal to 500 MeV. 

A great advantage of the EFT approach to nuclear forces is that it creates two- and many-body forces on an equal 
footing. Three-nucleon forces make their appearance at the third order in the chiral power counting. These leading- 
order contributions are: the long-range two-pion exchange graph; the medium-range one-pion exchange diagram; and 
the short-range contact term. 

In Ref. j!7j . density-dependent corrections to the in- medium NN interaction have been derived from the leading- 
order chiral 3NF. These are effective two-nucleon interactions that reflect the underlying three-nucleon forces and are 
therefore computationally very convenient, whereas realistic models of three-nucleon forces would be prohibitive. 

A total of six one-loop diagrams contribute at this order. Three are generated by the two-pion exchange graph of 
the chiral three-nucleon interaction and depend on the low-energy constants 01,3,4, which are fixed in the NN system 
[§]. We use C\ = —0.81 GeV -1 , c 3 = —3.2 GeV -1 , and c 4 = 5.4 GeV -1 . Two are generated by the one-pion exchange 
diagram and depend on the low-energy constant cu- Finally, the short-range component depends on the constant Ce- 
The constants cd and eg can be fixed by fitting properties of few-nucleon systems, such as the triton and 3 He binding 
energies [13]. We use cp = 5.0 and ce = 0.48. 

In pure neutron matter, the contributions proportional to the low-energy constants C4,C£>, and Ce vanish |17) . 
Analytical expressions for these corrections are provided in Ref. |17j in terms of the well-known non-relativistic two- 
body nuclear force operators. These can be conveniently incorporated in the usual NN partial wave formalism and 
the conventional BHF theory. 

IV. RESULTS 

In Fig. |3j we display the EoS of SNM obtained with the two approaches outlined above. The (blue) solid lines 
denote the results from the meson-exchange Bonn potential and the conventional BHF approximation ( "BnB BHF" ) 
or the relativistic DBHF ("BnB DBHF"), while the chiral calculations are shown by the (red) dashed lines, with and 
without 3NF, as denoted. 

As to be expected, both the conventional BHF calculation and the one using only the two-body chiral potential 
N 3 LO display excessive attraction and are unable to produce saturation up to very high density. The "Dirac effect" 
brings in a powerful saturation mechanism as seen by comparing the two (blue) solid curves. A very similar effect, 
both qualitatively and quantitatively, is generated by the effective chiral 3NF. 

In spite of the differences between the two approaches, the predictions are very close, not only in the value of the 
energy at saturation, but also the density dependence of the saturation mechanism. We notice, though, that the 
DBHF energies tend to grow faster at the higher densities. 

The same comparison is done in Fig. [4] for neutron matter, yielding similar conclusions. Here, too, the DBHF 
energies show a faster growth as density increases. 

In Fig. 5, we display the symmetry energy as obtained with the well-known parabolic approximation, e S y m — 
ejvM — esNM- (Notice that at this point we only consider the versions of the BnB- and the N3LO-based calculations 
which have realistic saturation properties, see Fig. 3.) 

The slope of the symmetry energy, usually defined in terms of the L parameter, (L = 3po( -^j m ) Po ), is a measure of 
the pressure gradient between neutron and symmetric matter. A variety of experiments are aimed at constraining this 
important quantity pQ, which correlates nearly linearly with the neutron skin of neutron- rich nuclei. At the present 
time there seems to be a consensus from various experiments that the acceptable range of values of the symmetry 
energy and its slope at saturation are centered around e sym and L equal to 32.5 MeV and 70 MeV, respectively pQ. 
The values of the L parameter obtained with "BnB+DBHF" and "N3LO+3NF" are 69 and 61 MeV, respectively. 
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FIG. 3: Energy /particle in SNM. The solid (blue) lines are obtained with Bonn B + DBHF (upper curve) and Bonn B + 
BHF (lower curve). The dashed (red) lines show the predictions by the N 3 LO potential with (upper curve) and without (lower 
curve) chiral three-nucleon forces. 

Notice that the typical uncertainty on the L parameter is as large as 20-25 MeV pQ. 

How the differences/similarities noted above impact neutron star bulk properties is examined in Fig. 6. (Since our 
focal point is a comparison between two theoretical approaches, rather than a detailed calculation of /3-stable matter 
composition, we consider only neutron stars made of pure neutron matter.) We also note that, in order to extend the 
chiral predictions to the high densities probed by compact stars, we fit a three-parameter function, e(p)=ap + ftp 1 , 
to the "N3LO+3NF" predictions of Fig. [4] and use this ansatz to obtain the high-density EoS needed for neutron star 



The comparison in Fig. 6 shows a larger star maximum mass for "BnB DBHF" , most likely the result of larger 
repulsion of that model at high density. Overall, the predicted star bulk properties are only moderately different in 
the two approaches. 



We have considered two very different methods to approach the study of nucleonic matter: one based on a meson- 
theoretic potential and the Dirac-Brueckner-Hartree-Fock approximation; the other based on a high-precision chiral 
NN potential and chiral effective three-nucleon forces at NNLO. The predictions we have considered include: the EoS 
of nuclear and neutron matter, the symmetry energy and its slope, and the mass-radius relation in a neutron star. 

From our results, we conclude that the DBHF method is an excellent phenomenology capable of incorporating 
important many-body effects that are crucial to nuclear saturation. 

In both approaches, the effective 3NF is generated by one nucleon interacting with the Fermi sea. That is, in both 
cases we have effective two-nucleon interactions that reflect the underlying three-nucleon forces. No matter if this 
interaction proceeds via relativistic meson exchange or via chiral EFT forces, the results are very similar. This is 
reassuring and confirms that the two ways of describing nuclear forces are complementary. 



calculations. 



V. 



CONCLUSIONS 
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FIG. 4: Same as Fig. 3, but for neutron matter. 
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FIG. 5: The symmetry energy vs. density for the models corresponding to the two upper curves in Fig. 3-4. 
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